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ABSTRACT 

We present observations of the Rossiter-McLaughlin effect for the transiting exoplanet 
systems WASP-1, WASP-24, WASP-38 and HAT-P-8, and deduce the orientations of 
the planetary orbits with respect to the host stars' rotation axes. The planets WASP-24b, 
WASP-38b and HAT-P-8b appear to move in prograde orbits and be well aligned, having 
sky-projected spin orbit angles consistent with zero: A = -4.7 ± 4.0°, A = IStgo and A 

= —9.7^7 70 , respectively. The host stars have Tcs < 6250 K and conform with the trend 
of cooler stars having low obliquities. WASP-38b is a massive planet on a moderately 
long period, eccentric orbit so may be expected to have a misaligned orbit given the high 
obliquities measured in similar systems. However, we find no evidence for a large spin- 
orbit angle. By contrast, WASP-lb joins the growing number of misaligned systems 
and has an almost polar orbit, A = — 79.0^4 30. It is neither very massive, eccentric 
nor orbiting a hot host star, and therefore does not share the properties of many other 
misaligned systems. 

Key words: stars: planetary systems - stars: individual: WASP-1, WASP-24, WASP- 
38, HAT-P-8 - techniques: radial velocities 
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scope, operated on the island of La Palma jointly by Denmark, Finland, 
Iceland, Norway, and Sweden, in the Spanish Observatorio del Roque 
de los Muchachos of the Institute de Astrofisica de Canarias; and the 
HARPS spectrograph mounted on the ESO 3.6m at the La Silla Obser- 
vatory in Chile under proposal 084.C-0185. 



1 INTRODUCTION 

The process of exoplanet migration has been a hotly de- 
bated topic since the first close-in planets were discovered 
jMavor & Ouelozll995l : rButler & Marcvll996l : lMarcv & Butleil 
Il996h . We are able to explore the mechanisms which move plan- 
ets inwards from large orbits through a statistical analysis of 
their dynamical properties such as orbital period, eccentricity 
and spin-orbit alignment. In particular, transiting planets allow 
us to measure the sky-projected angle (A) between the stellar 
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rotation axis and planetary orbit t iirough m easureme nt of the 
Rossiter-McLaughlin (RM) effect l lRossiteil fl924: Mc LaughlinI 
Tiie effect is caused by the planet sequentially passing 
over and blocking portions of the rotating stellar surface result- 
ing in a radial velocity (RV) shift which traces the trajectory of 
the planet across the stellar disc and allows the orbital obliquity 
to be estimated. 

The first measurements of A found the systems to be 
well aligned and suggested that the planets had lost or- 
bital angular momentum through interactions with the proto- 
planetary disc ( lLinetal.lll993 : iMurrav et al.ll 199^ . However, 
there are now a growing number of planets with highly mis- 
aligned and even retrograde orbits which now make up approx- 
imately one-third of the systems so far studied. Several theo- 
ries have been postulated to explain this. For example, it has 
been suggested that proto-planetary discs may not always be 
aligned with the stellar rotatio n axis as previously assu med 
teate. Lodato. & P ringle 20 1^; iLai, FoucMt. & Lidlioioh . al- 
though [^tsonetalj C 201 Obi) fo und no evidence fo r misal igned 
debris-discs in eig ht systems. Ipabrvckv & WinnI ( l2009h and 
IXriaud et al.l ( 1201 Ol ) suggest that another, more dynamically vi- 
olent process involving interactions with a third body (another 
plane t or star), causes the misaligned orbits (see Rasio & Ford 
Il99d : IWu & Murray! I2OO3I : iNagasawa et alj|2008ir Whether a 
single mechanism or a combination of several is at work re- 
mains to be tested by increasing the number of measured sys- 
tems. 

Misaligned orbits appear to be synonymous with ec- 
centricity; six of the eight eccentric systems with measured 
spin-orbit angles are reported to have obliquities significantly 
different from zero XO-3 (Hebrard et al. 2008), HP 80606 
I Moutou et al]l2009l:IPont et al., 2009. : .Gillon 200^:lwinn et al 



2009| ) . WASP-8 dOueloz et al.l|201(]|). WASP-14 jjohnson et al 



2OO9I), HAT-P-1 1 /Winn et al.' '2010c'), HAT-P-14 ("W innetal 



2010b) with HAT-P-2 (Loeillet et al. 2008) and HD 17156 
iNarita et al.ll2009 ) being the exceptions. Althoug h HD 17156 
does n ot show a large sky-projected misalignment. ISchlaufmanI 
( l2010h notes that the stellar rotation axi s may be tilted alon g 
the line-of-sight. Another trend, noted bv ljohnson et al.l ^20091) , 
is the correla tion between planet mass and misalignment. 
lHebrardetal] ( l2010l) suggests that there could be several popu- 
lations of planets, with those more massive than Jupiter under- 
going a diffe rent migration scen ario leading to th e hig h spin- 
orbit angles. IWinnetalJ bOlOah and ISchlaufmMil |2m0) find 
that misaligned orbits are more common in host stars with larger 
masses/higher effective temperatures. It is suggested that the 
tidal torques experienced by cooler stars, with deeper convec- 
tive zones, could cause their envelopes to quickly align with a 
planet's orbit, thereby erasing any initial misalignment. Excep- 
tions to this effect may be longer period and low mass planets 
which experience weaker tidal forces and therefore longer tidal 
timescales. 

We present spectros copic observat ions of the tran- 



siting planets WA SP- lb ( Cameron et al.l 12 007). WASP-24b 



dStreet et alJbOlQ) WAS P-38b dSarros etal..,2010) and HAT- 
P-8b jLatham et al.ll2009h obtained using the HARPS, SOPHIE 
and FIES spectrographs to determine their spin-orbit align- 
ments. All the host stars are of similar temperature, ~ 6100 K, 
but have diverse physical and dynamical properties which al- 
lows us investigate possible trends independent of T^a- WASP- 
I, WASP-24 and HAT-P-8 are non-eccentric, short period plan- 
ets (2.3-3.1 d) with masses similar to Jupiter (0.9-1.5 Mj), 



whereas WASP-38 is a massive planet on a moderately longer 
period orbit with a small but sig nificant eccentricity (mp = 2.7 
Mj, P = 6.9 d, e = 0.032+g °g^°. lBarros et eXXlOW ). 

In Section [2] we describe the general data analysis proce- 
dures and methods performed to measure the spin-orbit angles. 
The derived parameters for the four systems are presented in 
Section[3] and the implications of these are discussed in Section 

s 



2 DATA ANALYSIS 

2.1 Radial velocity extraction 

The data presented in this paper were obtained using three spec- 
trographs; HARPS, SOPHIE and FIES, the instrumental setup 
and data reduction for each is described in this section. Details 
of the individual observations and results of each system are 
discussed in Section[3] 



2.1.1 HARPS 

The HARPS instrument is a high resolution (R = 110,000) 
stabailed echelle spectrograph mounted at the La Silla 3.6 m 
ESO telescope. Observations were conducted in the OBJO 
mode without simultaneous Thorium- Argon (ThAr) calibration. 
The wavelength solution was calculated using a ThAr calibra- 
tion at the start of the nig ht, and HA RPS is stable within 1 
m s~^ across a night iRu pprecht et al . 2004), which is much 
lower than the photon noise on the data points. Spectra were ex- 
tracted and cross-correlated against a template of a G2V-type 
star using the HARP S Data Reduction Software (DRS), see 
iBaranne et al.' (1996), Pepe et al.! ( l2002h . iMavor et alj ( l2003h 
and ,Lovis & Pepe, C2007,) for more details. 



2.1.2 SOPHIE 

SOPHIE is a cross-dispersed, environmentally-stabilised 
echelle spectrograph (wavelength range 3872.4-6943.5A) 
designed for high-precision radial velocity measurements. The 
spectrograph was used in high efficiency mode (resolution R 
= 40,000) with the CCD in slow read-out mode to reduce the 
read-out noise. Two 3 arc-second diameter optical fibres were 
used, the first centred on the target and the second on the sky 
to simultaneously measure its background in case of contam- 
ination from scattered moonli ght. The spectra were reduced 
using the SOPHIE pipeline dPerruchot et alj l2008l) . Radial 
velocities were computed from a weighted cross correlation of 
each spectrum with a numerical mask of G2V spectral type. A 
Gaussian was fitted to the cross-correlation functions to obtain 
the radial velocity shift, and the uncertaint y was com puted 
using the empirical re lation given bv iouchveTSl lIoO^ and 
ICameron etall j2007h . 



2.1.3 FIES 

The FIES spectrograph is mounted on the 2.5-m Nordic Opti- 
cal Telescope on La Palma. FIES was used in medium resolu- 
tion mode (R = 46,000) with simultaneous ThAr calibration. 
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Spectra were extracted using the bespoke data reduction pack- 
age FIEStooQ An IDL cross-correlation routine was used to 
obtain the radial velocities by finding the maximum of the cross- 
correlation functions of 30 spectral orders and taking the mean. 

A template spectrum was constructed by shifting and co- 
adding the out-of-transit spectra, against which the individual 
spectra were cross-correlated to obtain the final velocities. This 
template was cross-correlated with a high signal-to-noise spec- 
trum of the Sun to obtain the absolute velocity to which the 
relative RVs were shifted. We estimated the RV uncertainty by 
a = RMS(t;) / v'iV, where v is the RV of the individual orders 
and A'^ is the number of orders. 

When fitting the data, we found that the FIES observations 
has a reduced Xred = X^^dof > 1 (dof = number of points - 
number of fitted parameters). This suggests that there is an ex- 
tra source of noise present which has not been accounted for 
in the internal errors. The cause may be that the ThAr calibra- 
tion does not travel th rough the same light path as the stellar 
light teuchhave et alJl201Q) . or that because the fibre lacks a 
scrambler there is n on-uniform illumination of the spectrograph 
jOuelozetalll 19991) . To account for these instrumental effects, 
15 m of uncorrelated noise needed to be added in quadra- 
ture to the internal error estimation to obtain a reduced value 
of unity. The photon errors are given in the data tables, whereas 
the rescaled uncertainties are shown in the figures. 
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2.2 Model fitting 

The RM effect and orbit were fitted simultaneously using 
all the available spectroscopic data. A Keplerian model was 
used for the or bit, and the analytical approach described in 
lOhta. Taruva. & Sutd dlOOSl) was used to calculate the RM ef- 
fect. We refer the reader to this paper for a detailed description 
of the model. The equations were modified to make them depen- 
dent on Rp/R, and a/i?, rather than Rp, i?, and a, to reflect 
the parameters derived from photometry and reduce the number 
of free parameters. 

We modelled the RM velocity anomaly as vb.m = 
—5vp/ (1 — 5) where 5 is flux blocked by planet ~ (Rp/R*)^ 
and Vp is sub-planet velocity, i.e. the velocity component of the 
rotating stellar surface blocked by the planet ^ Xp v sin i /R„, 
where Xp is the x co-ordinate of the position of the planet on the 
stellar surface, see Figure 5 of lohta. Taruva, & Sutd ( l2005h . 

In brief, the model comprises the following parameters: 
the orbital period, P; mid-transit time (in the UTC system). To; 
planetary to stellar radius ratio, Rp/ R-t \ scaled semi-major axis 
a/R,; orbital inclination, i; orbital eccentricity, e; longitude of 
periastron, uj; radial velocity semi-amplitude, K; systemic ve- 
locity of orbital dataset, 7; sky projected angle between the stel- 
lar rotation axis and orbital angular momentum vector. A; pro- 
jected stellar rotational velocity, vsini; and the stellar linear 
limb-darkening coefficient, u. 

Each dataset was allowed to have a different systemic ve- 
locity (7) to account for instrumental offsets. A linear limb dark- 
ening law was assumed, as the quadratic law alters the model by 
only a few m s~^ and so does not seem justified given the pre- 
cision of the RM data. As a test, u was left as a free parameter 
and no significant effect on A or u sin i was found. 



http://www.not.iac.es/instruments/fies/fiestool/FIEStool.html 



Figure 1. Top: Photometry of three transits of WASP-lb and the resid- 
uals from the best-fit model. 



Some parameters have been tightly constrained by previ- 
ous observations (e.g. P, i?p/i?t). We use this information in 
the form of a penalty function on the statistic; 

r 1 

2 X ^ ^i.obs ^z.calc 

f A - Aobs + aAobs X G(0,l) y 

V Aobs / 

where Hi.obs and Hi.caic are the ith observed and calculated ra- 
dial velocities and at is the corresponding observational error. 
A is a fitted parameter, Aobs is the parameter value determined 
from other observations and UAobs is the uncertainty in Aobs. 
The value G(0, 1) is a Gaussian randomly generated number of 
mean = and standard deviation = 1 . This allows the uncertainty 
in Aobs to be accounted for in the error budget. Depending on 
the situation, several penalty functions were used to constrain 
parameters in the fit and these are described in the individual 
analyses. If a parameter has asymmetric uncertainties, we have 
adopted the larger error. 

Best-fit parameters were obtained by minimi sing the 
statistic using the IDL-based MPFIT function ( iMarkwardtl 
l2009l) ; a least-squares minimisation technique using the 
Levenberg-Marquardt algorithm. The la best-fit parameter un- 
certainties were calculated using a Monte-Carlo method. We 
created synthetic data sets by adding a la Gaussian ran- 
dom variable to the data points. The free parameters were re- 
optimised for each simulated data-set to obtain the distribution 
of the parameter values. These distributions were not assumed 
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Figure 2. Left: Phase folded radial velocities of WASP-1 minus the systematic velocity (given in Table|5), over-plotted with the best-fit model, with 
the residuals shown below. The orbital observations were taken using FIES (open triangles) and SOPHIE (open circles) and the transit sequence using 
SOPHIE (filled circles). Right: Spectroscopic transit minus the orbital velocity, over-plotted with the best-fit model (A = — ,vsmi = 
5.77 it 0.35 km s~^, and residuals shown below. The dotted hne represents the RM effect of an aUgned orbit. 



to be Gaussian and the \a limits were found from the 15.85% 
and 84.15% bounds. 



3 ANALYSIS 



3.1 WASP-1 



A transit of WASP- lb was observed with the SOPHIE spec- 
trograph at the 1.93m telescope at Haute-Provence Observatory 
on the night of 2009 September 24. We acquired 21 spectra of 
WASP-1 covering the full transit (224 minutes, 12 observations) 
and a period of duration 156 minutes (9 observations) post tran- 
sit. No observations were possible prior to transit due the the 
low altitude of the target at this time. The exposure time was 
adjusted between 13 and 22 minutes in order to reach a constant 
signal-to-noise ratio (S/N) of 30 at 550 nm. Atmospheric condi- 
tions were stable and the exposure time did not require altering 
to obtain the required S/N. Moon illumination was 39% and at a 
distance of > 115° so this did not significantly affect the radial 
velocity determination. 

To fit the orbit, we used 7 S OPHIE observations o btained 
during the discovery of the planet jCameron et al.l2007l) and 16 
FIES spectra taken at various orbital phases between 2009 Jan- 
uary 1 and 2009 September 27. T he new SOPHIE and F IES 
observations are shown in Table [T] I wiieatlev et al.1 lioToh ob- 
served a secondary eclipse of WASP- 1 and found that the ec- 
centricity was consistent with zero, and in the models, e has 
been set to accordingly. The linear limb-darkening coefficient 
was chosen from the tables of lClared ( 12004) (ATLAS models) 
for the g' filter and fixed at ii = 0.73. 

A significant period of time has passed since the discov- 
ery of WASP-1 in 2006 and the orbital ephemeris may have 
drifted, leading to an inaccurate determination of the mid-transit 
time. To upd ate the ephemeris, we fitte d the z-band light curve 
prese nted in Charbonneau et d] j2007h . the full I-band transit 
from lShporer et al' l l2007h and a high quality light curve taken 
on 2010 August 26 by F. Harmuth as part of the TRESCA/ETD 



project ('Poddan v. Brat. & Peichj|2010 i). as shown in Figure[T] 
The TRESCA transit was observed using the 1.2-m telescope at 
Calar Alto observatory, Spain using an R-band filter. A Markov 
Ch ain Monte Carlo (MCM C) routine was used to fi t the data, 
see ICameron et al.l ( |2007h and IPollacco et alj j2008l) for more 
details. The non-linear limb darkening coefficients were chosen 
for the appropriate stellar temperature and photometric pass- 
band for each light curve. The parameters found for P, To, 
Rp/Rt, a/Rf and ip were used to constrain the fit in the form 
of penalty functions, and are shown in Equation|2l 

WASP- lb has a low impact parameter, b < 0.2, and in this 
regime v sin i and A are highly correlated. The shape of the RM 
signal is not strongly dependent on A, whereas the amplitude is 
dependent on both usini and A. To break this degeneracy we 
introduced a penalty function on v sin i using the value of 5.79 
± 0.35 km s ~^ found from spectros copic line broadening mea- 
surements bv lStempels et alJ ( l2007l) . The following statistic 
was adopted: 



+ 



vsini- 5.79 km s~^ + 0.35 km s~^ x G(0, 1) 
0.35 km s"^ 



P - 2.5199449 d + 0.0000013 d x G(0, 1] 
0.0000013 d 



+ 



+ 



To - 54461.86099 + 0.00024 x G{0, 1) 
0.00024 

Rp/R, - 0.10271 + 0.00059 x ^(0, 1) 
0.00059 

g/fl, - 5.64 + 0.12 X G(0, 1) \ ^ 

0.12 ) ^ 

ip- 88.07° + 1.15° xG(0,l] 
1.15° 



http://var2.astro.cz/EN/tresca/ 



+ 



(2) 
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Figure 3. Left: Phase-folded radial velocities of WASP-24 over-plotted with the best-fit model. The systematic velocity has been removed (see Table 
|6j and the residuals are shown below. The orbital points were taken with FIES (open triangles), CORALE (open squares) and HARPS (open squares) 
and the HARPS transit sequence is shown in filled circles. Right: The RM effect minus the orbital velocity is shown over-plotted with the best fit model 
and that of an aligned orbit (dotted line). The residuals are plotted below. 



The measured parameters and uncertainties are given in 
Table [5] Figure |2] shows the data over-plotted with the best-fit 



model, A 



-79.0^ 



and the model corresponding to an 



aligned orbit, A = 0° . It is obvious that the amplitude and sym- 
metry of the observations do not match with the aligned model. 
We conclude that WASP- 1 has a severely misaligned orbit with 
respect to the rotation axis of the host star. 

The value of v sin i is constrained by the penalty function 
and thus the uncertainty in A is small despite the low impact 
parameter. It must be noted, however, that deriving v sin i from 
spectroscopic line broadening is susceptible to systematic er- 
rors due to the uncertainties in broadening mechanisms, such as 
micro- and macro turbulence, which are difficult to quantify. By 
altering v sini by ± 1 km s~^ we find that the best-fit value of 
A changes by no more than 5°, so the interpretation of a highly 
obliquity is retained. The major factor leading the high precision 
in the determination of A is the system geometry. In a system 
with a low impact parameter and large misalignment, the form 
of the RM effect changes rapidly with A so it easy to differenti- 
ate between different angles. (jaudi__& Winn 12()07) provide an 
relationship to estimate the expected uncertainty in A. Substi- 
tuting ctrv = 8.5 m s~^, A'obs = 12 and the derived values of 
Rp/Rt, b, V sin i and A into their Equation 16, we obtain ax = 
5.8°. This is in good agreement with the value we obtain and 
the small difference may due to the assumption that v sin i is 
a free parameter whereas here we constrain it with the penalty 
function. 



FIES and 18 CORALIE out -of-tran sit observations taken during 
the discovery of the planet jStreet et al..201Q) . The HARPS ob- 
serv ations are given in Table|2] 

IStreet et aP (1201 Ol) reported that the orbit is not eccentric so 
e was fixed to 0. The l inear limb-da rkening coefficient was cho- 
sen from the tables of lClarell | |2004) (ATLAS models) for the V 
filter and fixed at = 0.66. Values of the photom etric param- 
eters P , To, Rp/R*, a/Rt and ip are taken from lStreet et al.l 
and are constrained through the fitting statistic be- 
low: 



E 



'^'i.obs '^i.calc 



+ 



P - 2.3412124 d -f 0.0000020 d x G(0, 1) 

0.0000020 d 
To - 55081.37941 + 0.00017 x G(0, 1) 

0.00017 

Rp/R, - 0.10040 + 0.00060 x G(0, 1) 
0.00060 



a/R, - 5.90 + 0.33 x G(0, 1) 
0.33 

ip - 83.64° +0.31° X G(0, 1) 
0.31° 



+ 



+ 



+ 



(3) 



3.2 WASP-24 

A transit of WASP- 24b was observed with the HARPS spec- 
trograph on the night of 2010 April 10. We acquired 51 spec- 
tra covering the full transit (161 minutes, 22 observations) and 
a period of duration 212 minutes (28 observations) distributed 
evenly before and after transit. The atmospheric conditions were 
stable and the seeing was 0.7". We also obtained 2 observations 
on the nights prior to and post transit to constrain the 7 velocity 
of the transit dataset. 

To fit the orbital parameters we obtained 8 HARPS obser- 
vations between 2010 March 26 and 2010 March 28 and used 10 



Figure[3]shows that WASP-24 has a very syrrmietrical RM 
effect, moving from red shift to blue shift. This implies that that 
the planet moves in a prograde, well-aligned orbit. A fit to the 
observations indicates an obliquity to be A = —4.7 ± 4.0°. All 
the fitted parameters are shown in Table (6) The fitted value of 
V smj = 7.32 ± 0.88km s"^ matches very well wit h that found 
from spectral line fitting, ■i;sini=7.0± 1.0 km s ~ ^ jStreet et"al] 
l2010h . WASP-24b has a high impact parameter (6 = 0.65) so a 
penalty function on v sin iwas not needed. We found that adding 
such a penalty fun ction constraining v sin i to the value from 
IStreet et^ l l201(i) had no significant effect on the derived pa- 
rameters. 
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Figure 4. Left: Phase folded radial velocities of WASP-38 over-plotted with the best-fit model. The systematic velocity has been removed (see Table 
[7) and the residuals are shown below. The orbital points were taken with SOPHIE (open circles) and CORALE (open triangles) and the FEES transit 
sequence is shown in filled circles. Right: The spectroscopic transit is shown minus the orbital velocity and over-plotted with the best fit model. The 
point represented by an open triangle was omitted from the analysis (see text). The dotted lines represent the models for A = 0°, 20°, 40°, 60° and 
80° (from top to bottom) showing the change in the shape of the RM effect with A. The residuals from the best-fit model are shown below. 



3.3 WASP-38 

The FIES spectrograph was employed to observe a transit 
of WASP-38 on 2010 June 8. We obtained 27 spectra dur- 
ing the night using an exposure time of 900s, giving 16 in- 
transit and 12 out-of-transit observations. The atmospheric con- 
ditions were stable and the seeing was 0.6". To fit the orbit 
we used the SOPH IE and CORALIE out-of-transit data from 
iBarros et al] bOlOh . The CORALIE radial velocity point at B JD 
= 2455404.6205 was removed, as it occurred during transit. To 
obtain Xred = 1, we required 15 and 8 m s^^ to be added in 
quadrature to the SOPHIE and CORALIE internal uncertainties, 
resp ectively. 

iBarros et al.l bOld) found that WASP-38 has a small but 
significantly eccentric orbit, so we fixed e = 0.0321 and u) = 
— 19°, since those parameters have little effect o n the fitted 
parame ters once the transit ephemeris is specified iWinn et ^ 
l2010bl). The photometric parameters found from light curve 
fitting tearros et alj|201ol) are constrained through the penalty 
function shown in Equation |4] We determined the linear limb- 
darkening coefficient from the tables of Claret (2004) (ATLAS 
models) for the V filter and fixed u = 0.64. 

The impact parameter of WASP-38 is low and has a rela- 
tively large uncertainty, 6 = 0.27to ij". As with WASP-1, we re- 
quire an independent constraint on v sin i to break the degener- 
acy between v sin i and A. The values of the photo metric param- 
eters a nd V sin i = 8.6 ± 0.4 km s~^ are taken from lBarros et alj 
( l2010h and used in the \^ statistic: 



Rp/R, ~ 0.0844 + 0.0011 x G(0, 1) 
0.0011 

a/R» - 12.15 + 0.19 x G(0, 1) 
0.19 



+ 0.67° X G(0, 1] 



+ 



0.67° 

Although the precision of the radial velocity measurements 
is lower than from the other spectrographs, we can see from the 
symmetry and characteristic red- then blue-shift of the RM ef- 
fect that the planet is generally aligned and not retrograde. A 
fit to the data yields A = IS^^go which is consistent with zero; 
the other fitted parameters are listed in Table |7] The large un- 
certainty in A is due in part to the small differences in the shape 
of the RM effect when A and h are small. Angles greater than 
approximately 60° produce a much more noticeable change in 
shape per degree than at low angles, as shown in Figure|4] Thus 
the data allow us to rule out very high misalignment angles, 
however further observations are needed to reveal whether a 
small misalignment exists. 

One observation lies more than three sigma from the best- 
fit model, see the open triangle in Figure |4] We investigated 
whether this could be due to an incorrect wavelength calibra- 
tion, but this does not appear to be the cause. However, seeing 
changes or guiding issues may be factors. The fit is improved by 
removing the point, so it was excluded in the analysis. The best- 
fit value of A obtained when the discrepant point is retained is 
— 42^*^30 which is not inconsistent with zero at the 1.3(T level. 



(4) 



= E 



+ 



vsini - 8.69 km s"^ + 0.40 km s"^ x G(0, 1) 

0.40 km s"^ 
P - 6.871815 d + 0.000045 d x G(0, 1) 
0.000045 d 

' To - 55335.92050 + 0.00074 x G(0, 1) 
0.00074 



3.4 HAT-P-8 

We observed a transit of HAT-P-8 with the FIES spectrograph 
on 2010 August 31. An exposure time of 900s was used to 
obtain 14 in-transit and 14 out-of-transit spectra. The seeing 
was 0.55" and weather conditions were good. Between the 
sixth and seventh in-transit observations, the telescope was re- 
pointed when the rotator reached the maximum limit, which re- 
quired ~9 minutes but did not affect the quality of the dataset. 
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We used the HIRES out-of-transit radial velocity points from 
iLatham et alJ ( 2009) to constrain the orbit and fixed the eccen- 
tricity of the orbit to zero. The linear limb-darkening coefficient 
was set to the value found from the tables of Claret ( 2004 ) (AT- 
LAS models) for the V filter, u = 0.69. 

Discovery light curves of HAT-P-8 were taken in 2007, so 
in order to improve the accuracy of th e ephemeris, we fitted the 
two full z-band light curve presented in Latha m et al ]||2009') and 
a high quality transit taken on 2010 August 28 by F. Harmuth 
(TRESCA/ETD), as shown in Figure [5] The TRESCA obser- 
vations were taken using the 1.2-m telescope at the Calar Alto 
Observatory in the R-band. As with WASP-1, an MCMC rou- 
tine was used to fit the data and the parameters found for P, To, 
Rp/Rt, a/Rt and ip are shown in Equation|6] 

An initial fit to the RM data found wsini ~ 16 km s^^ 
(and A ~ 8° ) to be significantly larger than the value of v sin i = 
11.5 ± 0.5 km s~^ reported by Latham et al. ( 2009). Other sim- 
ilarly fast rotating stars have also been found to show a discrep- 
ancy between the value of v sin i derived from sp ectroscopic 
line broadening and from RM measurements (e.g. IWinn et al] 
l2007l : iTriaud et al.ll2009l ). It is likely to be the result of the as- 
sumption that the line profile asymmetry caused by the planet 
blocking the rotating stellar surface can be modelled as a shift 
in the mean line position. For faster rotating stars, the asymme- 
try in the line profile is better resolved because of the broader 
line width, wh ich causes a larger apparent shift than expected, 
as seen here. Iffirano et alj fcoOft) have addressed this issue 
by mo difying the equations presented in lOhta. Taruva. & Sutd 
( 120051) to c ompensate fo r this effect. We implemented this so- 
lution as in lSimpson et a l. ( 2010b) and Bavli ss et aL C2010,) by 
calculating the RM radial velocity shift as: 



V-RM 



-S Vr, 



-1 3/2 



1 + 



2^2 + Cr2 



2/32 + a2 



(5) 



where /3 is the intrinsic line width and a — v sin i/a, where a 
is a scaling factor depending on limb darkening parameters (see 
Equation F6 of H09). We determined /3 as the quadrature sum 
of the macrotubulence (I'mac ), microturbulence (umic ) and in- 
strumental profile. Following IValenti & Fischej ( l2005h we used 
Vmac = 4.60 km s~^ and Wmic = 0.85 km s~^. The instrumental 
profile is given by c/R = 6.38 km s~^, where R is the instru- 
mental resolution 46,000. Thus we determined j3 = 7.9 km s^^. 
We calculated a = 1.31 for the limb darkening coefficients 
ui — 0.69 and U2 = 0. 

As with WASP- lb and WASP-38b, HAT-P-8b has a low 
impact parameter and therefo re we placed a penalt y function on 
the value of v sin i found by iLatham et al] (|2009|) to constrain 
the fit. The statistic used was: 



= E 



'^'i.obs "^i.calc 



+ 



?;sini - 11.5kms~i +0.5kms-i x G(0, 1) 
0.5 km s"^ 

P - 3.0763350 d + 0.0000040 d x ^(0, 1) 
0.0000040 d 



+ 



To - 54437.67582 + 0.00034 x G(0, 1) 
0.00034 

0.09110 + 0.00090 X 0(0,1) 
0.00090 



+ 



^ Rp I R» 



0.00 



0.15 



Lalhom el. ol (2009) z-bond 
_ 2007 Ocl 29 , 



Lolhom el, ol (2009) z-bond 
2007 Dec 02 



F. Hormuth (TRESCA) R-band 
2010 Aug 28 



' - V<^- '' '^''-'"' '' ""'^*«'-'' -' ''V , vV^ * . 



-0.04 



-0.02 



0.00 
Phose 



0.02 



0.04 



Figure 5. Photometry of three transits of HAT-P-8 and the residuals 
from the best-fit model. 



ajR, - 6.35 0.34 x (7(0,1) 

0.34 ' ^ 



ip - 87.26° + 1.00° X (7(0,1) 
1.00° 



(6) 



From the form of the RM effect, HAT-P-8 appears to have 
a prograde aligned orbit, as shown in Figure [6] A fit to the data 
yields A = — 9.7^y'^o , and the other fitted parameters are shown 
in Table m As with WASP-38b, the precision of FIES and the 
combination of low h and A limits the accuracy of the result, but 
allows us to rule out highly misaligned orbits. A transit of HAT- 
P-8b was independently observed with SOPHIE (Moutou et al., 
in prep.); it provides results similar to the ones presented here. 

Two observations show a large deviation from the best-fit 
model and again we could not determine an instrumental cause. 
We found that the fit was not significantly altered nor improved 
by removing the points so chose to retain them. The best-fit 
value of A having removed the two points (-18± 11°) is con- 
sistent with that found if the points are retained. 



4 CONCLUSIONS 

The spectroscopic transits of WASP- lb, WASP-24b, WASP-38b 
and HAT-P-8b have been observed using the HARPS, SOPHIE 
and FIES spectrographs. We modelled the RM effects and found 
the sky projected spin-orbit alignment angle of the systems to 
be A = -79.0+*1o, A = -4.7 ±4.0°, A = 15+^?° and A = 



-9.7+^;?° respectively. WASP-24b, WASP-38b and HAT-P-8b 
do not appear to be strongly misaligned and are consistent with 



© 0000 RAS, MNRAS 000, 000-000 



8 E. K. Simpson et al. 




Orbitol phose Orbitol phose 

Figure 6. Left: Phase-folded radial velocities of HAT-P-8 minus the systematic velocity (given in Table [D and over-plotted with the best-fit model 
with residuals shown below. The orbital observations were taken using HIRES (open circles) and the transit sequence using FIES (filled circles). 
Bottom: The spectroscopic transit, minus the orbital velocity, is shown over-plotted with the best-fit model with the residuals shown below. A dotted 
line represents the RM effect of an aligned orbit. 



zero within 2g. This suggests they could have undergone a rela- 
tively non-violent migration process which did not perturb them 
from the primordial alignment of the proto-planetary disc. Al- 
ternatively, tidal interactions may have forced the s tellar spin 
to align with the planetary orbit jWinn et ai]|2010j) . By con- 
trast, WASP- lb joins the approximately one-third of planets 
with misaligned orbits and points towards a dynamically vio- 
lent evolution such as gravitational scattering by other planets 
or the three-body Kozai mechanism. 

The uncertainties in A for HAT-P-8b and particularly 
WASP-38b are relatively large and allow for the possibility that 
the systems have a small misalignment. These larger uncertain- 
ties are due, in part, to the geometry of the systems. Both plan- 
ets have low impact parameters and in this regime the shape of 
the RM effect, which largely determines A, changes very subtly 
with A until angles of >60° are reached. Thus, many values of 
A can be fit by the data and the uncertainty is large. This ex- 
plains how we are able to obtain a much more precise value for 
WASP- lb which has a value of |A| which is much greater than 
zero. In addition, we are also limited by the instrumental noise 
on the FIES spectrograph and we recommend that further obser- 
vations of these systems be undertaken to refine the parameters. 

We modelled the RM effect using the equations of 
lOhta. Taruva. & Sutd JioO^ . However, it has been shown that 
this formulation does not estimate v sin i well for more rapidly 
rotating stars. We found this to be the case for HAT-P-8, which 
was measured to have v sin i = 11.5 ± 0.5 km s"^ from spec- 
troscopic line broadening but a best-fit value of 16 km s^^ from 
th e RM effect. For thi s case, we used the modified formulation 
of lHiranoet"aD(l2009l) to model the data, and placed a constraint 
on V sin i using the value from spectroscopic line broadening to 
break the degeneracy between A and v sin i at low impact pa- 
rameters. These modifications did not affect the conclusion that 
HAT-P-8b appears to be well aligned. 

The RM effect allows us to measure the spin-orbit angle in 
the plane of the sky. However, the line-of-sight alignment can be 
estimated by comparing the inclinations of the planetary orbit 
and stellar rotation axes ip and i, . Transiting planets naturally 
provide a tight constraint on ip given that the orbit must be close 
to edge-on in order for a eclipse to be observed. The stellar spin 



axis is more difficult to obtain and may be inferred by compar- 
ing the line-of-sight rotational velocity, v sini, to the true stel- 
lar rotation velocity v. Theoretical predictions of v are provided 
bv lSchlaufinai](l201(]h for many of the host stars of transiting 
planets. However, it can also be determined through the relation 
V — 2iiR,/Prot, where Prot, the stellar rotati on period, can 
be measured directly from sp ot modulation (e.g. lBaliunas et 
1 19971 : ISimpson etatlllOlOah or indirectly from Ca II H & K 
measurements ('Noves et al. 1 1984 IWatson et al] |20 1 Oah if this 



information is available. 

The four stars in this study have similar va l ues o f 
Tcff (see Tables |5] - [8]l and age (see ISchlaufmanI |2010|) , 
so we would expe ct them to have comparable values of v. 
ISchlaufmanI (2010) predicts this to be of the order of 7-9 
km s~^. In practice, we find a wide range of values of v sini; 
5.77 ± 0.35 km s-\ 7.32 ± 0.88 km s-\ 8.58 ± 0.39 km s"^ 
and 11.8 ±0.5 km s-\ This can be explained by supposing 
that the slower rotating stars are tilted so that we only measure a 
fraction of the true rotational velocity. WASP-1 is slower rotat- 
ing than p redicted (v sin i = 5. 8 km s^^, Vsi-m = 8.6 km s~^) and 
from this, ISchlaufmanI ( |2010|) deduced that the star and planet 
are significantly misaligned in the line-of-sight. This method 
yields four solutions for the inclination of the rotation axis of 
the star, u ~ 40°, 140°, 220° and 320°. 

The total alignment angle can be calculated through the 
following relation: 



cos tp — COS i, cos ip + sin i, sin ip cos A. 



(7) 



Substituting ip = 88.6° and A = -79°, we obtain V = 82°, 
84°, 96° and 98°. We therefore conclu de that WASP-lb has an 
almost polar orbit. ISchlaufmanI ilOldt did not find WASP-24b 
and HAT-P-8b to have highly misaligned orbits in the line-of- 
sight, and WASP-38b has a similar v sin i, so we do not expect 
it to be misaligned either. This agrees with the lack of evidence 
for strong misalignments in the plane of the sky from the RM 
effect in these systems. 

Previous authors have noted correlations between spin- 
orbit misalignment and eccentric orbits and massive planets. 
With e = 0.032 and rup = 2.7 M,j, WASP-38b may be expected 
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to be misaligned, however we find no evidence for high spin- 
orbit mgles^^ 

IWinn et alj \WlOa) noted that planets orbiting hot stars, 
defined as Tefr> 6250 K, tend to have misaligned orbits, 
whereas stars cooler than this host aligned planets. Low mass 
and long period planets are exceptions, as the timescale for tidal 
alignment is longer and may not have taken effect. All the plan- 
ets in this study have temperatures placing them in the cool 
star category and so WASP-24b, WASP-38b and HAT-P-8b ap- 
pear to conform with this trend. However, with its ~ 7 d orbit, 
WASP-38b has the potential to be an exception and to have re- 
tained its primordial spin-orbit angle. 

By contrast, the WASP-1 system, with T^ti = 6110 ± 45 
K, is unexpectedly misaligned. In this way, WASP-lb is similar 
to CoRoT-lb, which also orbits a cool star, Tcft- 5950 ± 150 
K tearge et alJlzOO S). and has a highly misaligned, almost po- 
lar orbit, A = 77 ± 11° (Pont et al. 2010) yet is not eccentric. 
Neither planets are low mass nor long period so do not meet the 
criteria to be exceptions from aligned orbits around cool stars. 
To bring the stars into the hot-star regime, the temperatures of 
CoRoT-1 and WASP-1 would have to be several a hotter than 
their best fit value, respectively. Whilst this could be the case, 
another possibility is that the distinction between cool and hot 
stars is not discrete and there may be a transition region between 
the two. 

Interestinglv. fcameron et all ((2007) reported a nearby stel- 
lar companion to WASP-1, and if associated, would have a pro- 
jected orbital separation of ~1800 AU. In this regime the Kozai 
mechanism can take effect and it is thought to significantly per- 
turb planetary orbits. Although beyond the scope of this paper, 
it would be very informative to investigate whether the migra- 
tion of WASP-lb could have been caused by interactions with a 
binary companion, if it is indeed associated. 

The rapid increase in the number of measured spin-orbit 
angles, and the discovery of a large population of misaligned 
systems, has revolutionised our understanding of planetary mi- 
gration. It is no longer believed that planet-disc interactions can 
fully explain the observed distribution of angles and other mech- 
anisms must play a significant role. In this paper we have pre- 
sented four new measurements of spin-orbit angles, and we find 
that WASP-lb has a near polar orbit, despite predictions that it 
would be aligned as it orbits a cool star. By contrast, WASP-38b 
is a massive planet on a moderately long period, eccentric orbit 
so has a high likelihood of misalignment but it does not show 
evidence for a large obliquity. In order to fully understand the 
processes which are at work, further measurements are needed 
to all ow theories to be compar ed with a strong observational 
basis. iMorton & Johnsoi] i20ld) estimate that a total of 40-100 
measurements will be required to be confident in the mecha- 
nisms at work, and with the steady flow of new transiting planet 
discoveries, uncovering the mystery of migration may soon be 
within our reach. 
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Table 1. Radial velocities of WASP-1 measured with SOPHIE during 
transit and FIES at other orbital phases. 



BJD 


RV 


EiTor 


-2 400 000 


(kms-i) 


(kms-i) 


Planetary transit: 






55099.33177 


-13.2101 


0.0083 


55099.34602 


-13.2153 


0.0086 


55099.35950 


-13.2128 


0.0087 


55099.37319 


-13.2040 


0.0086 


55099.38674 


-13.2146 


0.0086 


55099.40012 


-13.2218 


0.0085 


55099.41520 


-13.2341 


0.0085 


55099.42765 


-13.2228 


0.0085 


55099.43928 


-13.2454 


0.0086 


55099.45150 


-13.2457 


0.0087 


55099.46349 


-13.2614 


0.0087 


jjuyy.4/j3z 


-IJ.zo/ 1 


U.UUoo 


55099.48773 


-13.2527 


0.0092 


55099.50089 


-13.2662 


0.0088 


55099.51269 


-13.2596 


0.0083 


55099.52912 


-13.2711 


0.0072 


55099.54277 




0.0081 


55099.55493 


-13.2635 


0.0086 


55099.56709 


-13.2843 


0.0087 


55099.58288 


-13.2739 


0.0088 


55099.59579 


-13.2883 


0.0087 


HJD 


RV 


EiTor 


-2 400 000 


(kms-i) 


(kms-i) 


Other orbital phases: 






54834.3378 


-13.3780 


0.0105 


54835.4419 


-13.6246 


0.0107 


54836.4229 


-13.4224 


0.0096 


55025.7036 


-13.4250 


0.0191 


55041.6483 


-13.5199 


0.0249 




1 ^ ^ifit^ 

- i J . JZUO 


0154 


55086.5259 


-13.4365 


0.0393 


55087.4126 


-13.6088 


0.0223 


55096.5906 


-13.4449 


0.0215 


55097.6051 


-13.6009 


0.0181 


55098.5769 


-13.4039 


0.0223 


55099.6120 


-13.5478 


0.0125 


55100.5339 


-13.5002 


0.0102 


55119.4279 


-13.4457 


0.0131 


55119.6984 


-13.5493 


0.0187 
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Table 5. Derived system parameters and uncertainties for WASP-1. Tlie effective temperature is taken from lStempels et aiH2007l) . 

Parameter (units) Symbol Value 

Free parameters: 
Projected alignment angle (° ) 
RV semi-amplitude (km s~^) 

Systemic velocity of SOPHIE transit dataset (km s~^) 
Systemic velocity of SOPHIE orbital dataset (km s^^) 
Systemic velocity of FIES orbital dataset (km s^^) 

Parameters controlled by priors: 
Projected stellar rotation velocity (km s~^) 
Period (days) 

Transit epoch (HID - 2 400 000) 
Planet/Star radius ratio 
Scaled semi-major axis 
Orbital inclination (°) 



Fixed parameters: 






Eccentricity 


e 





Limb darkening 


u 


0.73 


Effective temperature (K) 




61 10 ±45 



Table 6. Derived system parameters and uncertainties for WASP-24. The effective temperature is taken from lStreet etal] j2010l) . 



Parameter (units) 


Symbol 


Value 


Free parameters: 






Projected alignment angle (°) 


A 


-4.7±4.0 


Projected stellar rotation velocity (km s~^) 


V sin j 


7.0 ± 0.9 


RV semi-amplitude (km s~^) 


K 


0.1482 ± 0.0025 


Systemic velocity of HARPS transit dataset (km s^^) 


71 


-17.7871 ± 0.0017 


Systemic velocity of HARPS orbital dataset (km s~^) 


72 


-17.7905 ± 0.0019 


Systemic velocity of FIES orbital dataset (km s^^) 


73 


-17.9050 ± 0.0073 


Systemic velocity of CORALIE orbital dataset (km s~^) 


74 


-17.8031 ± 0.0040 


Parameters controlled by priors: 






Period (days) 


p 


2.3412127 ± 0.0000020 


Transit epoch (HID - 2 400 000) 


To 


55081.37941 ±0.00017 


Planet/Star radius ratio 




0.10040 ± 0.00060 


Scaled semi-major axis 


a/Rt 


5.98± 0.15 


Orbital inclination (°) 


i 


83.64 ± 0.29 


Fixed parameters: 






Eccentricity 


e 





Limb darkening 


u 


0.66 


Effective temperature (K) 


Tcff 


6075 ± 100 



A -79.0 

K 0.1094 ±0.00564 

71 -13.2365 ± 0.0027 

72 -13.5016 ± 0.0065 
7a -13.4990 ± 0.0057 



vsini 5. 77 ±0.35 

P 2.5199449 ±0.0000013 

To 54461.86099 ± 0.00024 

Rp/R* 0.10271 ± 0.00058 

a/Rt 5.64±0.13 
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Table 7. Derived system parameters and uncertainties for WASP-38. The effective temperature is taken from lBarros et alj )2010h . 



Parameter (units) 


Symbol 


Value 


Projected alignment angle (°) 


A 


-43 


RV semi-amplitude (km s~^) 


K 


0.2538 ± 0.0035 


Systemic velocity of FIES transit dataset (km s~^) 


71 


y.O'lUI _o.oo57 


Systemic velocity of SOPHIE orbital dataset (km s~^) 


72 


-9.7181 ±0.0063 


Systemic velocity of CORALIE orbital dataset (km s^^) 


73 


-9.7951 ±0.0027 


Parameters controlled by priors: 






Projected stellar rotation velocity (km s~^) 


V sini 


8.58 ± 0.39 


Period (days) 


P 


6.871814 ± 0.000045 


Transit epoch (HJD - 2 400 000) 


To 


55335.92044 ± 0.00074 


Planet/Star radius ratio 


Hp/ i?* 


0.0844 ± 0.0011 


vj^alCU. aCllll lllttlUl tlAlf* 


a/ Rf 




Orbital inclination (°) 


i 


88.83 toJs 


Fixed parameters: 






Eccentricity 


e 


0.032 


Longitude of periastron 


U! 


-19.0 


Limb darkening 


U 


0.64 


Effective temperature (K) 




6150 ± 80 



Table 8. Derived system parameters and uncertainties for HAT-P-8. The effective temperature is taken from lLatham et al] )2009l) 



Parameter (units) 



Symbol 



Value 



Free parameters: 

Projected alignment angle (°) A 

RV semi-amplitude (km s~^) K 

Systemic velocity of FIES transit dataset (km s~^) 71 

Systemic velocity of HIRES orbital dataset (km s~^) 72 



-9.7 



+9.0 



-7.7 

0.1580 ± 0.0041 
-22.3650 ± 0.0040 
-0.0937 ± 0.0031 



Parameters controlled by priors: 
Projected stellar rotation velocity (km s^^) 
Period (days) 

Transit epoch (HJD - 2 400 000) 
Planet/Star radius ratio 
Scaled semi-major axis 
Orbital inclination (°) 



V sin I 
P 
To 
Rpl R* 

a/Rt, 



11.8± 0.5 

3.0763370 ± 0.0000036 

54437.67587 ± 0.00034 

0.09135 ± 0.00089 

6 12+0'20 
-0.21 

fO.75 



87.80^ 



-0.77 



Fixed parameters: 
Eccentricity 
Limb darkening 





0.69 



Effective temperature (K) 



T„t 



6200 ± ! 
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Table 2. Radial velocities of WASP-24 measured with HARPS during 
and outside transit. 



BJD 


RV 


Error 


-2 400 000 


(km s ^1 


(km s^i) 


Planetary transit! 






55296.6514 


-17.7358 


0.0059 


55296.6651 


-17.7463 


0.0126 


55296.6701 


-17.7492 


0.0118 


55296.6750 


-17.7407 


0.0113 


55296.6801 


-17.7504 


0.0112 


55296.6850 


-17.7733 


0.0119 


55296.6900 


-17.7487 


0.0119 


55296.6950 


-17.7436 


0.0119 


55296.7000 


-17.7492 


0.0124 


55296.7051 


-17.7629 


0.0119 


55296.7100 


-17.7551 


0.0123 


55296.7150 


-17.7619 


0.0123 


55296.7201 


-17.7668 


0.0121 


55296.7250 


-17.7521 


0.0140 


55296.7300 


-17.7364 


0.0131 


55296.7350 


-17.7379 


0.0126 


55296.7400 


-17.7434 


0.0123 


55296.7450 


-17.7512 


0.0124 


55296.7501 


-17.7645 


0.0118 


55296.7550 


-17.7403 


0.0116 


55296.7600 


-17.7633 


0.0119 


55296.7650 


-17.7531 


0.0120 


55296.7700 


-17.8006 


0.0122 


55296.7751 


-17.7885 


0.0122 


55296.7800 


-17.7813 


0.0115 


55296.7851 


-17.7965 


0.0111 


55296.7901 


-17.8128 


0.0109 


55296.7950 


-17.8160 


0.0107 


55296.8000 


-17.8399 


0.0098 


55296.8049 


-17.8265 


0.0116 


55296.8100 


-17.8495 


0.0121 


55296.8150 


-17.8366 


0.0121 


55296.8199 


-17.8141 


0.0139 


55296.8251 


-17.7952 


0.0130 


55296.8300 


-17.8240 


0.0129 


55296.8350 


-17.7943 


0.0122 


55296.8401 


-17.8418 


0.0118 


55296.8451 


-17.8177 


0.0112 


55296.8500 


-17.8012 


0.0109 


55296.8550 


-17.8154 


0.0108 


55296.8600 


-17.8018 


0.0110 


55296.8650 


-17.8293 


0.0110 


55296.8701 


-17.8289 


0.0109 


55296.8750 


-17.8200 


0.0110 


55296.8801 


-17.8355 


0.0105 


55296.8851 


-17.8203 


0.0108 


55296.8901 


-17.8294 


0.0107 


55296.8951 


-17.8398 


0.0100 


55296.9000 


-17.8602 


0.0111 


55296.9051 


-17.8438 


0.0116 


55296.9101 


-17.8287 


0.0116 


Other orbital phases: 






55281.6933 


-17.7475 


0.0052 


55281 9080 


-17.6632 


0.0045 


55283.6945 


-17.8572 


0.0056 


55283.9139 


-17.7825 


0.0045 


55297.7382 


-17.8696 


0.0042 


55297.9076 


-17.8057 


0.0043 


55282.6762 


-17.7674 


0.0125 


55282.6835 


-17.7461 


0.0118 


55282.6911 


-17.7431 


0.0118 


55282.6988 


-17.7443 


0.0141 


55295.6485 


-17.7747 


0.0056 


55295.9092 


-17.6802 


0.0039 



Table 3. Radial velocities of WASP-38 measured with FIES during tran- 
sit. The starred point (*) was omitted from the analysis 



HJD 


RV 


Error 


-2 400 000 


(kms-i) 


(kms-l) 


55356.4021 


-9.8108 


0.0120 


55356.4135 


-9.8175 


0.0084 


55356.4248 


-9.8113 


0.0091 


55356.4362 


-9.8320 


0.0057 


55356.4476 


-9.7805 


0.0114 


55356.4590 


-9.7907 


0.0069 


55356.4704 


-9.8042 


0.0081 


55356.4817 


-9.7962 


0.0074 


55356.4931 


-9.8042 


0.0103 


55356.5045 


-9.8028 


0.0089 


55356.5159 


-9.8149 


0.0121 


55356.5273 


-9.8188 


0.0010 


55356.5386 


-9.8072 


0.0099 


55356.5500* 


-9.7912 


0.0098 


55356.5614 


-9.8347 


0.0074 


55356.5728 


-9.8802 


0.0075 


55356.5842 


-9.8762 


0.0106 




-9 8841 


0087 


55356.6070 


-9.9084 


0.0089 


55356.6183 


-9.8734 


0.0098 


55356.6297 


-9.8284 


0.0081 


55356.6411 


-9.8365 


0.0076 


55356.6525 


-9.8599 


0.0102 


55356.6639 


-9.8602 


0.0096 


55356.6752 


-9.8593 


0.0083 


55356.6866 


-9.8661 


0.0093 


55356.6980 


-9.8855 


0.0122 
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Table 4. Radial velocities of HAT-P-8 measured with FIES during tran- 
sit. 



HJD 


RV 
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-2 400 000 


(km s ) 
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